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Abstract 

We report on a measurement of the proton structure function F2 in the range 
3.5 X 10"^ < X < 4 X 10-3 and 1.5 GeV^ < < 15 Qev2 at the ep colhder 
HERA operating at a centre-of-mass energy of ^/s = 300 GeV. The rise of F2 with 
decreasing x observed in the previous HERA measurements persists in this lower x 
and range. The evolution of F2, even at the lowest and x measured, is 
consistent with perturbative QCD. 
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1 Introduction 



The measurement of deep inelastic scattering (DIS), ep eX, at HERA has shown a 
rapid rise of the proton structure function F2{x,Q'^) with decreasing a; for a; < 10~^ ^ ||. 
The corresponding increase of the virtual photon-proton cross section cr^o/ with the centre- 
of-mass energy W is much stronger than that of cr7o? fo^^ I'sal photons [Q, ^]. The slower 
increase of the real photoproduction cross section is consistent with the energy behaviour of 
hadron-hadron total cross sections. In perturbative QCD the rise of F2 at low x is ascribed 
to an increase in the sea quark density |^ and the significance thereof is discussed 
extensively in the literature [|^. One of the important questions is how far perturbative QCD 
retains its validity as one probes large parton densities. Until now the HERA measurements 
of F2 have covered values above 4.5 GeV^. It is of great interest to extend the F2 
measurement to lower values and study the low x behaviour in the transition region 
between photoproduction and DIS. 

In this paper, F2 measurements, from e^p DIS data, at x > 3.5 x 10^^ and at values 
as low as 1.5 GeV^ are reported. Access to such low values is achieved in two different 
ways. For the first analysis (SVX) HERA is operated with the interaction point shifted in 
the proton direction in order to improve the acceptance for small positron scattering angles. 
The second analysis (ISR) uses DIS events from the nominal interaction point with initial 
state photon radiation where the radiated photon is detected. These events effectively have a 
lower initial positron beam energy and thus for a fixed acceptance of the positron scattering 
angle events with smaller values of can be reconstructed. The reconstruction of events 
having positrons scattered at small angles is improved compared to the 1993 analysis by 
the addition of the small angle rear tracking detector (SRTD). 

2 Experimental Conditions 
2.1 HERA Running Conditions 

The data were taken with the ZEUS detector at HERA in 1994. HERA operated with 153 
colliding bunches of 820 GeV protons and 27.5 GeV positrons. Additional unpaired positron 
and proton bunches circulated, which are used to determine beam related background. The 
root mean square of the proton bunch length was approximately 20 cm while the positron 
bunch length was negligible in comparison, leading to an interaction length having a root 
mean square of 10 cm. For the SVX analysis, the mean interaction vertex was moved from 
Z = to Z = Q7 cm Q Approximately 5% of the proton current was contained in satellite 
bunches, which were shifted by 4.8 ns with respect to the primary bunch crossing time, 
resulting in a fraction of the ep interactions occurring on average +72 cm upstream of the 
primary position. 

^The ZEUS coordinate system is defined as right handed with the Z axis pointing in the proton beam 
direction, and the X axis horizontal, pointing towards the centre of HERA. The origin is at the nominal 
interaction point. 
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The SVX analysis is based on an integrated luminosity of 58 nb~^ collected while HERA 
operated with the shifted interaction point. The ISR analysis uses 2.5 pb~^ of data taken 
with HERA operating with the interaction point at the nominal position, Z = 0. 



2.2 The ZEUS Detector 

A description of the ZEUS detector can be found in 0, §. The primary components 
used in these analyses are the uranium-scintillator calorimeter (CAL) and the tracking 
detectors. The calorimeter covers 99.7% of the total solid angle and is subdivided into 
electromagnetic (EMC) and hadronic (HAC) sections with typical cell sizes of 5 x 20 cm^ 
(10 X 20 cm^ in the rear calorimeter (RCAL), i.e. in the positron beam direction) and 

20 X 20 cm^ respectively. The calorimeter has an energy resolution of a/E = 18%/ ^E{GeV) 



for electrons and a / E = 35%/ \ E{GeV) for hadrons, as measured in test beams. The timing 



resolution of a calorimeter cell is better than at = 1.5/y _E'(GeV) © 0.5 ns (© denotes 
addition in quadrature). 

The tracking system consists of a vertex detector (VXD) UTO) and a central tracking 
chamber (CTD) |TT[ enclosed in a 1.43 T solenoidal magnetic field. The interaction vertex 



is measured with a resolution along (transverse to) the beam direction of 0.4 (0.1) cm. 

The position of positrons scattered close to the positron beam direction is determined 
by the SRTD which is attached to the front face of the RCAL. The SRTD consists of two 
planes of scintillator strips, 1 cm wide and 0.5 cm thick, arranged in orthogonal directions 
and read out via optical fibres and photomultiplier tubes. It covers the region of 68 x 68 cm^ 
in X and Y and is positioned aX Z = —148 cm. A hole of 20 x 20 cm^ at the centre of 
the RCAL and SRTD accommodates the beampipe. The SRTD is able to clearly resolve 
single minimum ionising particles (mip) and has a position resolution of 0.3 cm. The time 
resolution is better than 2 ns for a minimum ionising particle. 

The luminosity is measured via the positron-proton bremsstrahlung process, ep e'yp, 
using a lead-scintillator calorimeter (LUMI) |jl2| which accepts photons at angles < 0.5 mrad 
with respect to the beam axis. The LUMI photon calorimeter is also used to measure 
the energy and position of photons from initial state radiation in DIS events. It is po- 
sitioned at Z = — 107 m and, under test beam conditions, has an energy resolution of 
a/E = 18%/ JE (GeV). In its operating position it is shielded from synchrotron radia- 



tion by a carbon-lead filter and has an energy resolution of a/E = 26.5%/-y/£'(GeV), as 
determined from bremsstrahlung data. The position resolution is 0.2 cm. In addition, an 
electromagnetic calorimeter positioned at Z = —35 m is used for tagging scattered positrons 
at small angles. 



2.3 Triggering 

Events were filtered online by a three level trigger system At the first level, DIS events 
were selected by requiring a logical AND between two conditions based on energy deposits 
in the calorimeter. The first condition was the presence of an isolated electromagnetic 
energy deposit. The EMC energy deposit was required to be greater than 2.5 GeV. The 
corresponding HAC energy deposit was required to be either less than 0.95 GeV or no more 
than a third of the EMC energy deposit. The threshold values have been chosen to give 
>99% efficiency for positrons with energy greater than 5 GeV as determined by Monte 
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Carlo studies. Further details of the algorithm can be found in |13|. The second condition 



required that the EMC section have an energy deposit greater than 3.75 GeV. Background 
from protons interacting outside the detector was rejected using the time measurement of 
the energy deposits from downstream veto counters and the SRTD. 

At the second level trigger (SLT), background was further reduced using the measured 
times of energy deposits and the summed energies from the calorimeter. The events were 
accepted if 

SsLT = T.E^il - cosOi) > 24 GeV - 2E^ (1) 

i 

where Ei and di are the energies and polar angles (with respect to the primary vertex 
position) of calorimeter cells, and is the energy deposit measured in the LUMI photon 
calorimeter. For perfect detector resolution and acceptance, 6slt is twice the positron 
beam energy (55 GeV) for DIS events while for photoproduction events, where the scattered 
positron escapes down the beampipe, 6slt peaks at much lower values. 

The full event information was available at the third level trigger (TLT). Tighter timing 
cuts as well as algorithms to remove beam-halo muons and cosmic muons were applied. The 
quantity Stlt was determined in the same manner as for 6slt- The events were required 
to have 

6tlt > 25 GeV - 2E^. (2) 

Finally, DIS events were accepted if a scattered positron candidate of energy greater than 
4 GeV was found. 



3 Monte Carlo Simulation 



Monte Carlo (MC) event simulation is used to correct for detector acceptance and smearing 



effects. The detector simulation is based on the GEANT programme [|14| and incorporates 
our understanding of the detector, the trigger and test beam results. Neutral current DIS 
events are simulated to 0{as) using the HERACLES programme which includes first 
order electroweak radiative corrections. The hadronic final state is simulated using the 



colour-dipole model including boson gluon fusion CDMBGF |T^ as implemented in ARI- 
ADNE [|r7| for the QCD cascade and JETSET ^ for the hadronisation. The ARIADNE 
model currently provides the best description of the observed DIS nondiffractive hadronic 
11 



final state 



20 



Diffractive events with a large rapidity gap as observed in the data 
are simulated within ARIADNE by assuming that the struck quark belongs to a colourless 
state having only a small fraction of the proton's momentum. The parameters of the model 
are adjusted to be consistent with recent ZEUS measurements [^. The MRSA |^ parton 
density parameterisations, modified at low as described in [^], are used. These param- 
eterisations provide an adequate description of previous ZEUS and HI measurements ||. 

The shape of the vertex distribution used in the simulation is taken from nondiffractive 
photoproduction events. For such events the vertex reconstruction efficiency is found to be 
high and only weakly dependent on the Z position of the interaction. 

For the shifted vertex analysis, a sample of events corresponding to an integrated lumi- 
nosity of ~ 100 nb~^ was generated with > 0.5 GeV^. 

The main source of background in the data sample for the SVX analysis is due to 
photoproduction leading to the detection of a fake scattered positron. Minimum bias pho- 
toproduction events are simulated using PYTHIA with cross sections according to the 
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ALLM parameterisation [25|. Photoproduction events corresponding to an integrated lumi- 
nosity of 110 nb~^ were generated with a photon-proton centre-of-mass energy W > 190 GeV. 
Events with smaller W values do not contribute to the photoproduction background. 

For the ISR analysis, a DIS sample corresponding to an integrated luminosity of 2.7 pb~^, 
with an initial state photon energy above 3 GeV, was generated with > 0.5 GeV^. 



4 Kinematic Reconstruction 

In deep inelastic scattering, e(/c) -|- p{P) e{k') + X, the proton structure functions are 
expressed in terms of the negative of the four-momentum transfer squared, Q^, and Bjorken 
X. In the absence of QED radiation, 

Q'- = -q' = -(k-k'f, (3) 

where k and P are the four-momenta of the incoming particles and k' is the four-momentum 
of the scattered lepton. The square of the centre-of-mass energy is denoted by s. The 
fractional energy transferred to the proton in its rest frame is y = Q^/{sx). 

In the SVX analysis, and x are reconstructed from the measured energy, E'^, and 
scattering angle, 9^ of the positron, 

Ql = AE,E'^cos\ej2), (5) 

^ E^Eicos\9j2) 
E,{E,-E'^s\n\ej2)y 

where Ee and Ep are the positron and proton beam energies and the scattered positron 
angle is measured with respect to the positive Z direction. 

The ISR sample is selected by requiring that a photon with energy E^ be detected in 
the LUMI photon calorimeter. The variables and x are determined using equations |^ 
and ^ replacing Ef. with E^. — E^, treating the virtual positron as a real positron, which is 
a good approximation for the ISR analysis. 



4.1 Vertex Determination 

The vertex coordinates are determined from tracks reconstructed with the CTD and VXD. 
The Z coordinate is determined on an event-by-event basis. Since the transverse sizes of 
the beams are smaller than the resolutions for the X and Y coordinates of the vertex, the 
beam positions averaged over the entire data sets are used. For events which do not have 
a tracking vertex, the Z coordinate is set to the primary position of the interaction point. 
For events with a tracking vertex, the resolution of Z is ±0.4 cm over the entire Z range 
considered in these analyses. At low y the current jet is produced at small forward angles, 
resulting in a reduced probability for vertex reconstruction. The vertex distributions for 
the two analyses, including events arising from the interaction of the satellite bunch, are 
found to be well reproduced by the MC simulation. 
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4.2 Positron Identification and Efficiency 



The positron identification algorithm is based on a neural network using information from 
the CAL and is described elsewhere The efficiency of finding the scattered positron is 



sensitive to details of the shower evolution, in particular to energy loss in material between 
the interaction point and the calorimeter. The efficiency was measured using elastic QED 
Compton events ep e'yp where the positron and photon are detected in the CAL by 
exploiting the presence of two electromagnetic objects and the overconstrained kinematics. 
DIS data were also used to study the positron finding efficiency by comparing different 
positron finding algorithms. A correction based on the QED Compton study is used in 
the analyses. The size of the correction to the efficiency obtained from the standard ZEUS 
detector simulation programme is 16% at 8 GeV decreasing to 0% at 18 GeV. The efficiency 
of the identification algorithm when the scattered positron has an energy of 10 GeV is 55%, 
rising to 100% above energies of 18 GeV. The uncertainty in the efficiency is accounted for 
in the systematic errors (see Sect. |5.5|). 



4.3 Positron Position Measurement 

The impact point of the scattered positron at the calorimeter is measured using the SRTD 
for the SVX data and part of the ISR data. The position resolution of the SRTD is deter- 
mined using positrons hitting the calorimeter at the boundary of cells where the position 
can be determined by the calorimeter with a resolution of less than 0.1 cm. The measured 
SRTD position resolution of 0.3 cm is well reproduced by the MC simulation. 

For events with a reconstructed vertex the scattering angle of the positron, determined 
from the event vertex and the impact position at the face of the SRTD, has a resolution of 
1.7 mrad, while for events without a tracking vertex the measured scattering angle has a 
resolution of 4 mrad. For the ISR analysis, 7% of the events are found outside the SRTD. 
For these events, the impact point is determined by the RCAL resulting in an angular 
resolution of 7 mrad. 



4.4 Positron Energy Measurement 

In the analyses, x and are determined using the corrected positron energy E'^ as described 
below. For background rejection and rejection of events with hard initial state radiation, 
we define, in addition, the quantity 6: 

6 = 6h + Se; 5, = ^E,(l-cos^,); 4 = ^^(1 - cos^J, (7) 

h 

where Eh is the energy deposited in the calorimeter cell h and the angle, 6h, is determined 
from the event vertex and the cell centre. The sum excludes the calorimeter cells belonging 
to the scattered positron. 6e is calculated using the corrected positron energy. 

In the range of the present analyses the scattered positron traverses typically two 
radiation lengths of passive material before reaching the calorimeter, thereby losing a sig- 
nificant amount of energy. The correlation between the energy lost in the inactive material 
in front of the calorimeter and the energy deposited in the SRTD is used to correct the 
calorimeter energy measurement. 

This correction can be determined from a data sample where the scattered positron 
energy is known from kinematic constraints. This is the case for DIS events with x = Ee/ Ep, 
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for which the scattered positron energy is equal to the incident positron beam energy. 
Events close to this condition, called kinematic peak (KP) events, can be selected by the 
requirement 6h < 0.06 ■ E^. The correction for lower energies can be obtained using QED 
Compton events. The energies of the positron and the photon can be predicted precisely 
from the measurement of their scattering angles assuming the transverse momentum of the 
scattered proton to be small. Events from elastic DIS p° production, ep — >■ epp"; p° 7r"''7r~ 
provide a way to check the SRTD energy correction. Here, the angle of the positron is 
determined from the event vertex and the position measurement at the face of the SRTD 
while the momenta of the vr^ and vr" are measured with the CTD. 

The SRTD energy correction is determined using the KP and QED Compton events. A 
clear correlation between the energy measured in the calorimeter and the energy deposited 
in the SRTD for KP events is observed (see Fig. [^a). Figure |I]b compares the corrected 
positron energies of the data with the MC simulation for KP events, where it can be seen 
that the peak and width of the data are well reproduced by the MC simulation. The 
deviation of the mean corrected energy from the prediction of the KP and QED Compton 
events is less than 2% as shown in Fig. |l|c for positron energies between 9 and 27.5 GeV. 
The point at the highest energy {E'^ = 27.5 GeV) is obtained from the KP events while the 
lower energy points are determined from the elastic QED Compton events. For the DIS p° 
data, the deviation of the mean corrected positron energy from the value expected from the 
p° TT+TT^ measurement is less than 1% (see Fig. |I|c); these data have not been used to 
determine the SRTD energy correction. The resolution of the corrected positron energy is 
shown in Fig. |l|d and can be described hj cr/E = 26% / E (GeV) . 

5 Analysis of the Shifted Vertex Data 

In the SVX analysis, F2 is measured using data from an integrated luminosity of 58 nb~^ 
collected with the vertex shifted by +67 cm in Z. The vertex shift extends the acceptance of 
positrons to smaller scattering angles and hence events with lower can be reconstructed. 

5.1 Event Selection 

The following cuts are used to select the DIS events for the SVX analysis. 

• The scattered positron energy as obtained from the calorimeter and corrected by the 
energy measured in the SRTD is required to be greater than 10 GeV. This ensures 
a high efficiency of finding the scattered positron and removes high y events which 
suffer from large photoproduction background. 

• The impact position of the scattered positron is required to be at least 3 cm from the 
inner edge of the rear calorimeter. This cut ensures that the electromagnetic shower of 
the positron is well contained in the calorimeter. The impact position of the scattered 
positron is required to be within the acceptance of the SRTD and at least 1 cm away 
from the edges of the SRTD. 

• The value of 6 for the event is required to be within 35 < 5 < 65 GeV in order 
to reduce photoproduction and beam-gas related background. This requirement also 
removes events with hard initial state radiation, thus reducing the radiative correc- 
tions. The hadronic energy measurement affects the event selection via this cut when 
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(5e < 35 GeV. The simulation of 6h in the MC is sensitive to the hadronisation scheme 
and the details of the hadronic energy loss in the inactive material. The mean 6h was 
compared to 6^ as a function of 6e and the MC simulation was found to reproduce the 
data adequately. 

• For events with a tracking vertex, the Z coordinate of the vertex is required to be 
within 25 cm < Z < 200 cm, the acceptance being extended to larger Z values to 
accommodate the events from the satellite bunch. If no vertex was reconstructed, the 
Z coordinate is set to the primary interaction point. 

A total of 13210 events pass the selection cuts. 

The distributions of the reconstructed x and for the selected events are shown in 
Fig. ^ The dashed line (Z=0, 1993 cuts) gives the limit of the acceptance for the scattered 
positron corresponding to the 1993 F2 measurement using the ZEUS detector. The accep- 
tance limit in this analysis is shown as the solid line (Z=67 cm, 1994 cuts). The region 
between the dashed and solid lines shows the increased acceptance obtained for the SVX 
data and by using the SRTD. 

5.2 Resolution of the Kinematic Variables and Bin Selection 

The selected events are binned as shown in Fig. 0. The sizes of the bins are determined by 
the resolution in x and Q^. Within the selected bins, the resolution in is found to be 
9-12% and the systematic shifts are less than 5% as determined from MC simulation. At 
the lowest values of x, the resolution in x is 20% increasing to 85% for larger x values where 
larger bins are chosen to ensure a high purity (defined as the number of events measured 
in the bin which originated from the bin divided by the number of events measured in the 
bin). The average reconstructed values of x are typically shifted by 6%, with shifts of up 
to 28% occurring at the lowest y values. The bias in the bins results from the steepness of 
the distribution and the radiative events. 

The purity of the events in the bins is typically about 45% and is greater than 30% for 
all bins. The acceptance (defined as the number of measured events originating from the 
bin divided by the number of events generated in the bin) is typically 70% except for the 
lowest bin where it is about 22%, and for the bin with 1.9 GeV^ < < 2.7 GeV^ 
and 0.9 xlO"^ < a; < 2.0 xlO"^ where it is 27%. In these bins the purity is 50% and 36% 
respectively. 

5.3 Background Estimation 

The background from beam-gas related interactions is estimated from events due to the 
unpaired bunches of positrons and protons. Events triggered from empty bunch crossings 
are used to estimate the background from cosmic ray showers. The surviving background 
events are subtracted statistically. The (x, Q^) bin with the highest background has 2% 
background. In the other bins no events survive and the background is estimated to be less 
than 1%. 

The largest contributor to background is photoproduction, where an energy deposit in 
the calorimeter is misidentified as a scattered positron. This background is estimated in 
two ways. The first method uses the photoproduction MC simulation. The same selection 
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procedure as for DIS events is applied and the accepted photoproduction events are sta- 
tistically subtracted for each (x, Q^) bin. The second method uses the 6 distribution to 
estimate the background. The two estimates are found to agree within statistical errors. 
The photoproduction background is found to be significant only in bins corresponding to 
the highest y values where it amounts to 2%. 



5.4 Determination of Fo 



In the range of the present analysis, the effect from Z° exchange is negligible and the 
double differential cross section for single virtual photon exchange in DIS is given by 



2(1 + 



1 + R 



F2(x,g2) 1 + 5,,(x,Q2) 



(8) 



where R is related to the longitudinal structure function, Fl, hy R = Fl/{F2 — Fl). The 
correction to the Born cross section from radiative effects is given by 5r- In the analysis, 
hard radiation coUinear with the final state positron is not resolved. Furthermore, the 
event selection criteria (see above) remove events with hard radiation collinear with the 
beam positron. Thus, in the kinematic range covered in this analysis, the effective radiative 
correction is typically 10% or less and is, to a sufficient approximation, independent of the 
structure function F2. 

An iterative procedure is used to extract the structure function F2. In the first step 
a bin-by-bin correction obtained from the MC simulation using the parton distributions 



given by MRSAp2| is applied to the data. The result for F2 from this first iteration is used 



for a QCD fit using DGLAP equations in next-to-leading order (QCD NLO fit) very 
similar to that described in The evolution uses massless quarks of three flavours in the 
proton, and the charm quark coefficient functions from [2^ to ensure a smooth crossing of 
the charm threshold; the charm contribution to F2 is calculated in leading order only. The 
NMC data for > 4 GeV^ are used to constrain the fit at high x. The value of R was 
taken using the QCD prescription 
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and the parton distributions from the QCD NLO 
fit. The effect of the R correction on the F2 values is significant only in the highest y bins 
where it is as much as 7%. 

The resulting QCD NLO fit parameters are used to reweight the MC events, and the 
procedure is repeated, leading to a new estimate of F2. The procedure is repeated until 
the F2 values from two consecutive iterations change by less than 0.5%. The final result 
is reached in three iterations. It should be noted that the QCD NLO fit is used here only 
as a parameterisation to obtain a stable acceptance correction and not to perform a QCD 
analysis. 

The statistical errors of the F2 values are calculated from the number of events measured 
in the bins and the statistical error on the acceptance calculation from the MC simulation. 
Since we are using bin-by-bin corrections, the correlations of statistical errors between the 
F2 measurements enter only via the finite statistics of the MC sample. The correlations are 
small given the relatively large MC sample used in this analysis. A correlation between the 
F2 values of neighbouring bins is present due to the acceptance and smearing corrections. 
The sensitivity of the measured F2 to this effect has been checked by using the F2 obtained 
from the MRSA ||2^ parton density parameterisations, modified at low as described 
in for the acceptance correction calculations. The change in the measured F2 values is 
found to be small compared to the statistical errors. 
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The measured distributions of and x are shown in Figs. ^ and 03 respectively. The 
data are shown as the sohd circles and the MC simulation results, normalised to the lu- 
minosity of the data and reweighted by the QCD NLO fit, are shown as histograms. The 
measured distributions of the positron energy and angle are shown in Figs. ^ and ^d. 
All events which pass the selection criteria described in Sect. with a reconstructed 
> 1 GeV^ are shown, including those events which fall outside the bins used for the 
analysis. These distributions have not been corrected for background. There is adequate 
agreement between data and simulation for both the shape and normalisation of the dis- 
tributions. The number of events, and values of R and F2 at x and values specified, 
which are chosen to be convenient for comparisons with other F2 measurements, are given 
in Table |I[ The functional form of the QCD NLO fit was used to calculate F2 at these x 
and values. Using alternate parameterisations for this correction has a negligible effect 
on the quoted values of F2. 



5.5 Systematic Uncertainties 

The systematic uncertainty of the F2 values is determined by changing the selection cut or 
analysis procedure and taking the difference between the measured value of F2 and the new 
value. Positive and negative differences are separately added in quadrature to obtain the 
total systematic error. For each change, the photoproduction background estimated from 
the MC simulation is first subtracted. 

The systematic uncertainties are grouped in the following categories: 

Positron energy and finding efficiency 

• The parameters of the neural network positron finder are varied resulting in changes 
of F2 by typically less than 4%. In addition, an alternate positron finder, which has 
been used in previous analyses [|^, is used as a check of the neural network finder. 
Consistent results are obtained in the region where the efficiencies of both finders are 
reasonably high. 

• In the MC simulation the reconstructed positron energy is increased and decreased 
by a linear function (2% at 5 GeV and 1% at 27.5 GeV). The magnitude of the shift 
represents our present understanding of the energy scale. These shifts have a 5% effect 
on F2 for the low x bins and 10% for the highest x bins. 

• The positron identification efficiency is varied within the errors of its determination 
based on the QED Compton study. The effect on F2 is negligible except in the lowest 
X bins where it is as much as 3%. This also checks the effects of a possible energy 
dependence of the trigger efficiency. 



Positron scattering angle 

• Changing the fiducial cut of the positron position at the inner edge of the calorimeter 
from 3 cm to 4 cm has a 3% effect on F2 for the lowest x and bins. 

• The reconstructed Z vertex position in the MC simulation is shifted by ±0.4 cm in 
Z. The effect on F2 is less than 1.5% for all the bins. 
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• The positron position reconstructed by the SRTD is shifted by ±0.2 cm in X and 
±0.15 cm in Y. The size of the shifts are estimated from the uncertainty of the 
position of the SRTD relative to the beam position in X and Y. These shifts have a 
7% effect on F2 for the lowest x bins and less than 4% for the other bins. 

Additional systematic uncertainties 

• The photoproduction background estimate is changed by +100% and —50% resulting 
in a 2% effect on F2 for the lowest x bins and a negligible effect for the higher x bins. 

• In the calculation of the acceptance, the fraction of the total cross section arising from 
diffractive scattering is increased by 60%. The effect on F2 is less than 5%. 

• The efficiency for reconstructing a vertex is determined by MC simulation and is found 
to decrease from 85% at y = 0.7 to 40% at y = 0.03. For events at low y the tracks 
in the CTD are at low angles as the current jet is produced in the forward direction, 
resulting in a reduced probability of reconstructing the event vertex. In the data, 
70% of the events have a tracking vertex compared to 73% in the MC simulation. 
The vertex of the events without a tracking vertex is set to the nominal shifted vertex 
position. The reconstructed values are larger than the true values if the events 
without a tracking vertex originated from the satellite bunch. The effect on F2 of the 
satellite bunch was studied by MC simulation and it was found to be largest at high 
X values where it is 8% decreasing to about 1% for the low x points. 

As an additional check, the Z coordinate of the vertex is fixed at the primary inter- 
action point for all events in the data and in the MC simulation in order to check 
the sensitivity to the vertex distribution. The change in F2 is largest for the high 
values at large x where it is about 15% of F2. For the lowest x and bins the effect 
on F2 is about 10%, while for the other bins it is less than 5%. The change in F2 from 
this systematic check is not included when determining the systematic error. 

• The 6 cut is lowered from 35 GeV to 32 GeV and raised to 38 GeV. The largest effect 
on F2 of about 3% occurs at the low x and low bins, and it decreases with increasing 

and x. Changing the 6 cut checks the photoproduction background estimate, the 
simulation of the QED radiative effects and the simulation of the hadronic energy 
measurement. 

The total systematic error on F2 of this analysis is 4 to 14% to be compared with the 
statistical error of 4 to 9%, which includes the statistical error from the MC simulation 
(see Table |l]). In addition to the above errors, there is an overall normalisation uncertainty 
of 3% due to the uncertainty in the first level trigger efficiency and to the uncertainty on 
the determination of the luminosity. The given errors of the F2 data do not include this 
normalisation uncertainty. 

6 Analysis of the Radiative Events 

The ISR analysis is based on data, corresponding to an integrated luminosity of 2.5 pb~^, 
taken in 1994 with the interaction point at Z = 0. F2 is measured using events with QED 
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initial state radiation. The emission of a hard photon from the initial state positron results 
in a scattering at a reduced centre-of-mass energy, and thus for a given scattering angle 
events with lower are accepted. 

6.1 Event Selection and Bin Selection 

The sample of DIS events with initial state radiation is selected in a manner identical to 
the SVX analysis, with the following exceptions: 

• An energy deposit in the LUMI photon calorimeter is demanded to ensure collinear 
radiation. The energy is required to be between 6 and 18 GeV in order to reduce 
background and ensure sufficient resolution of the reconstructed kinematic variables. 
Events having an energy greater than 3 GeV measured in the electron calorimeter of 
the luminosity detector are rejected in order to reduce the accidental overlap of DIS 
and photoproduction with those from the bremsstrahlung process ep ejp. 

• The corrected energy of the scattered positron detected in the RCAL is required to be 
greater than 8 GeV. For positrons within the fiducial volume of the SRTD, the energy 
is corrected as in the SVX analysis. In areas not covered by the SRTD a position 
dependent description of the inactive material in front of the GAL is determined using 
the KP events and is used to correct the measured energy in the GAL. A description of 
the method can be found in |T|. A correction for the positron identification efficiency 
in the MG simulation is also made as in the SVX analysis. 

• The fiducial cut for the positron position is the same as that described in the SVX 
analysis for the inner edge of the SRTD. If the positron position is outside the fiducial 
volume of the SRTD, the position is reconstructed using the calorimeter. 

• The 6 cut is replaced by a cut on 6' = 6 + 2 ■ E^. 

A sample of 10726 events is selected. The resolutions of x and are similar to those 
given in Sect. |5.2| and depend only weakly on the photon energy E^. The bins in x and 
are chosen in a similar way to those used in the SVX analysis. 

6.2 Background Estimation 

The main source of background is the accidental overlap of DIS or photoproduction events 
with bremsstrahlung events, ep e'jp. This background is estimated with events selected 
with the above criteria, removing the requirement of a tagged photon and the S' cut. To 
these events a photon is added with an energy determined from a random sampling of 
the measured spectrum of photon energies for bremsstrahlung events. The resulting 6' 
spectrum is normalised to the 6' distribution of the data above 70 GeV where no DIS 
events are possible, taking energy and momentum conservation and the resolution of the 
detector into consideration, to obtain an estimate of the background. It is subtracted from 
the data for each {x, Q"^) bin and is below 10% of the number of events observed in the bin 
except in the lowest x bins where it is up to 24% in one bin. 

The beam related background is determined and subtracted statistically in the same 
manner as described in Sect. ^. It is below 5% for the three lowest bins except for the 
bin at the lowest x where it is 7%. It is negligible in higher bins. The background from 



11 



photoproduction with initial state radiation and cosmic rays is negligible. Event losses due 
to bremsstrahlung overlap with a DIS event having initial state radiation are also negligible. 



6.3 Determination of F2 

F2 is determined as described in Sect. p.4| . The MC sample is first reweighted using the 
QCD NLO fit to the F2 from the SVX analysis. F2 is then determined using the reweighted 
MC for the acceptance correction. A second QCD NLO fit to the measured ISR F2 values 
(excluding the SVX F2 values) is performed and the final F2 is determined. The measured 
distributions of and x are shown in Figs. |a and respectively. The data are seen to 
be in adequate agreement with the MC simulation after the background is added. 

The photon energy spectrum of the selected events without the cut on is shown in 
Fig. ^ and is well described by the expectation from the DIS MC simulation and back- 
ground estimates. Figure Ijd shows the difference between the energy of the radiated photon 
measured in the LUMI photon calorimeter and E^^^, the photon energy determined, using 
energy and momentum conservation, from the measurements of the energies and the angles 
of the scattered positron and hadrons in the CAL. The agreement between the MC simu- 
lation and the data is reasonable. The distribution is centred at zero indicating agreement 
in the energy scales of the LUMI photon calorimeter and the CAL. 



6.4 Systematic Uncertainties 

The checks used to estimate the systematic uncertainty are similar to those presented in 
Sect. except for the checks related to the displacement of the vertex. To estimate the 



uncertainties due to the photon tagging the following additional checks are made: 

• The lower cut on the photon energy is raised from 6 to 7 GeV. The effect on F2 is below 
2.5% for the bins at low y and around 6% for the bins at high y and is compatible 
with statistical fluctuations. 

• The resolution of the photon energy measurement in the MC simulation is degraded 
from 26.5% ■ ^E (GeV) to 28.5% ■ ^E{GeV). The effect on F2 is around 1.5%. 

• The energy calibration of the LUMI photon calorimeter in the MC simulation is 
changed linearly, within the estimated uncertainties, by 0.4% at the positron beam 
energy to 3% at 5 GeV. The effect on F2 is around 7%. 

• The acceptance of the LUMI photon calorimeter for ISR events obtained from the 
MC simulation is about 30%. An independent determination of the acceptance using 
bremsstrahlung data leads to an effect of about 4% on F2. 

The beam divergence is determined with the LUMI photon calorimeter using the 
bremsstrahlung events. If the beam divergence is increased by 15% in the acceptance 
calculation, the effect on F2 is, on average, 1.7%. 

• The event selection is repeated without the requirement that the LUMI electron energy 
be less than 3 GeV. The effect is below 3% for the bins at low y and is up to 10% in 
two bins at high y. 
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The total systematic uncertainty of this analysis is 14 to 27% to be compared to the 
statistical error of 7 to 14%. The effect of additional QED radiative corrections not included 
in the HERACLES MC programme is small relative to the measurement errors in the 



kinematic region considered here [BI 



7 Results and Discussion 

The measured F2 values from the SVX data (ZEUS SVX 1994) are listed in Table |I], and 
those from the ISR analysis (ZEUS ISR 1994) are listed in Table The data are presented 
in Fig. ^ together with previous ZEUS measurements (ZEUS 1993). The three data sets 
are in good agreement. The present analyses significantly increase the measured kinematic 
region. Using our data alone, the rise of F2 as x decreases is observed down to = 3 GeV^. 
Including the fixed target data (the F2 parameterisations shown as curves in Fig. ^ provide 
a good description of the fixed target data at higher x) the rise is seen to persist down to 
Q2 = 1.5 GeV^. 

In Fig. ||the measured F2 values are compared to the prediction of GRV(94) |^ which is 
based on perturbative QCD using the DGLAP evolution equations. The GRV(94) parton 
distributions have a very low starting scale, Qq, for the DGLAP evolution equation of 
0.34 GeV^, where the gluon and sea distributions are assumed to have valence-like spectra. 
The steep rise in F2 at low x is generated dynamically by the evolution in Q^. The GRV(94) 
predictions for F2 are in agreement with the data showing that perturbative QCD can 
describe the data down to values of 1.5 GeV^ at the low x values of this measurement. 



The predictions of Donnachie and Landshoff based on Regge phenomenology are also 



shown and are seen to be ruled out for >2 GeV^ and disfavoured for = 1.5 GeV^. 

The DIS cross section can be expressed as the product of the flux of virtual photons 
and the total cross section a^ol for the scattering of virtual photons on protons [0. al^t 
is defined in terms of the cross section for the absorption of transverse and longitudinal 
photons, (Jt and ai respectively, by 

atr^aT{x,Q^) + aL{x,Q^). (9) 

The expression for F2 in terms of (Tt and 0"^ is 

F,(x, = ^l^g^^;^[^T(x, Q') + a,{x, Q% (10) 

where rrip is the mass of the proton. The separation into the photon flux and cross section 
can be interpreted in a way similar to the interaction of real particles provided that the 
lifetime of the virtual photon is large compared to the interaction time, or x ^ 1/ {2mpRp) 
where Rp ^ 4 GeV~^ is the proton radius |3^. At small x the expression can be written 
in terms of the total virtual photon-proton centre-of-mass energy W (where W"^ = + 
Q2(i/x- 1)) to give 

The measured F2 data are converted to the total virtual photon-proton cross section and 
shown in Fig. |^ along with low energy data and real photoproduction cross section mea- 
surements. The DIS data are seen to rise steeply as a function of W'^ between the fixed 



atJ{W\Q^)^^F2{x,Q'). (11) 
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target and the HERA energy range, even at values as low as = 2.0 GeV^ in contrast 
to the cross section for real photons which exhibits only a slow rise from the fixed target 
data to the HERA data g i- 

In summary, the proton structure function F2{x, Q^) has been measured in DIS in a new 
kinematic range down to = 1.5 GeV^ and x = 3.5 x 10~^. The GRV(94) predictions, 
which are based on perturbative QCD, are found to be consistent with the data. The 
Regge predictions of Donnachie and Landshoff based on the soft pomeron are ruled out for 

> 2 GeV^ and disfavoured for = 1.5 GeV^. For the centre-of-mass energies between 
the fixed target regime and the HERA energy range, the total virtual photon-proton cross 
section aigf is found to rise steeply with the centre-of-mass energy for Q"^ as low as 2.0 GeV^. 
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Table 1: The measured F2{x,Q'^) from the SVX analysis. The bin boundaries and values 
of X and at which F2 is determined are listed. The numbers of events before background 

subtraction as well as the estimated photoproduction and beam-related background (in the 
column labeled "No. BG events") for each bin are given. The values of R, which are 
used to determine F2 from the differential cross sections, are shown (see text). An overall 
normalisation error of 3% is not included. 
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Table 2: The measured F2{x,Q'^) from the ISR analysis. The bin boundaries and values, of 
X and at which F2 is determined are listed. The numbers of events before background 
subtraction as well as the estimated accidental event overlap and beam-related background 
(in the column labeled "No. BG events") for each bin are given. The value of R obtained 
from the NLO Q CD fit is tabulated. An overall normalisation error of 3% is not included. 
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Figure 1: a) The correlation between the energy measured in the calorimeter and the energy 
in the SRTD in units of mips (mean energy deposited by one minimum ionising particle) for 
the kinematic peak (KP) events, b) The distribution of the corrected positron energy for the 
data, shown as the points, and the MC simulation, shown as the histogram, for KP events. 
The arrows in a) and b) indicate the positron beam energy, c) Measured fractional deviation 
between the mean corrected calorimeter energy and the predicted energy as a function of 
the corrected energy, d) The measured energy resolution of the calorimeter as a function of 
energy. The curve corresponds to 26% -^J E{GeV), the resolution used in the MC simulation. 
For c) and d), the data are results using QED Compton (squares), DIS elastic p° (triangles) 
and KP (dots) events. See text for details. 
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Figure 2: The x-Q"^ distribution for events passing the selection criteria from the 1994 SVX 
analysis. The extension in the accepted region compared to the 1993 analyses is shown 
between the dashed line labeled "Z=0, 1993 cuts" and the solid line labeled "Z=67 cm, 1994 
cuts (shifted vertex)". See text for more details. 
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ZEUS 1994 




Figure 3: a) The reconstructed Ql of the SVX event sample, b) The reconstructed Xe 
distribution of the SVX event sample, c) The spectrum of the scattered positron energy, d) 
The distribution of the positron scattering angles. In the figures the data (dots) are compared 
with the MC simulation (histograms). All events with a reconstructed Ql > 1 GeV^ which 
pass the selection criteria described in Sect. |5. j| are shown. The background has not been 
subtracted from the data. The MC distributions have been reweighted using the final F2 
parameterisation from the QCD NLO fit to the ZEUS data and normalised to the luminosity 
of the data. 
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ZEUS 1994 




Figure 4: a) The reconstructed Ql distribution of the ISR sample, b) The reconstructed 
Xe distribution of the ISR sample, c) The spectrum of the photon energy measured in the 
LU MI photon calorimeter without the cut on the photon energy, d) The difference between 
the photon energy measured in the LUMI photon calorimeter and that determined from the 
CAL. In the figures the data (dots), background estimate (hatched histogram) and sum of 
the background and DIS MC (solid histogram) are shown. All events with a reconstructed 
Ql > 1 GeV^ which pass the selection criteria described in Sect. 6.1 are shown. The MC 
distributions have been reweighted using the final F2 parameterisation from the QCD NLO 
fit to the ZEUS data and normalised to the luminosity of the data. 
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ZEUS 1994 




Figure 5: The measured F2 from the SVX analysis (solid dots), the ISR analysis (solid tri- 
angles) and the 1993 results (open squares) compared with the expectations from GRV(94) 

(solid line) and Donnachie and Landshoff (DL) (dashed line). Overall normalisation uncer- 
tainties of 3% for the 1994 results and 3.5% for the 1993 points are not shown. The inner 
error bars represent the statistical errors while the outer error bars represent the systematic 
errors added in quadrature to the statistical errors. 
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Figure 6: The total virtual photon-proton cross section versus W"^ for different values. 
The cross section values obtained from the F2 values described in this paper and the F2 values 
from the 1993 data are shown in addition to data from previous low energy experiments 
/p^. The region to the right of the dashed line correspond to x < l/{2mpRp). Also shown 
is the W"^ behaviour of the measured cross section for real photoproduction together with the 
prediction of Donnachie and Landshoff jS^/ (solid line). 
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